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ABSTRACT 


Electrostatic  plugging  of  multidipole  cusps  is  investigated. 


The  plugging  is  produced  by  placing  positively  biased  electrodes  on 


both  sides  of  all  magnetic  cusps.  Experimental  observations  show 


the  plasma  potential  follows  the  electrode  potential  and  the  plasma 


electron  density  increases  by  a  factor  of  about  two  as  the  electrode 


bias  increases  frcm  zero.  The  addition  of  a  negatively  biased  grid 


at  one  end  of  the  chamber  does  not  change  the  response  of  the  plasma 


electron  density  to  the  electrode  bias.  A  model,  which  considers 


plasma  losses  along  the  magnetic  cusps  as  well  as  by  diffusion  trans¬ 


verse  to  the  field  lines,  is  shown  to  qualitatively  agree  with  the 


data.  The  elevated  interior  plasma  potential  permits  the  identifi¬ 


cation  and  investigation  of  secondary  electron  production  at  the 


walls.  Data  is  presented  which  shows  that  secondary  electron  emission 


rate  at  the  chamber  walls  can  be  comparable  to  plasma  production  rate. 


A  measurement  of  the  relative  ion  densities  of  Hydrogen  is  attempted. 


but  the  results  are  inconclusive  as  to  the  effect  improved  confine¬ 


ment  has  on  the  various  ion  species  densities. 
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1 


I .  INTRODUCTION 

The  use  of  fields  to  confine  laboratory  plasma  has  been 

studied  for  many  years.  The  use  of  magnetic  fields,'*'  to  increase 

particle  path  lengths,  is  limited  by  losses  along  the  field  lines. 

For  example,  the  use  of  cusps  and  surface  magnetic  fields  has  shown 

that  the  loss  area  for  ions  is  bounded  by  Ur^L  (Ref.  3)  and  2r.  L 

He  i  c 

(Ref.  4)  where  r^  =  r^  ,  r^  is  the  ion  gyroradius,  rg  is  the 

electron  gyroradius,  and  Lc  is  the  length  of  the  magnetic  cusps. 

2 

Attempts  using  electrostatic  fields,  to  keep  either  ions  or  electrons 

from  the  chamber  walls,  have  failed  because  of  heating  of  the  grid 

wires  or  capacitor  plates,  used  to  produce  the  fields,  from  plasma 

bombardment  and  also  because  of  arcing  problems. 

To  improve  plasma  confinement,  various  authors  have  tried  a 

5 

combination  of  electrostatic  and  magnetic  confinement.  Dolan 
describes  this  technique  as  "...  magnetic  shielding  of  the  grid  of  an 
electrostatic  plasma  confinement  device  or  electrostatic  plugging  of 
a  magnetic  confinement  device."  Most  of  this  early  work  has  involved 
ring  spindle  cusp  devices,  with  electrodes  biased  negative  (with 
respect  to  the  plasma) .  The  present  work  is  the  first  attempt  at 
using  electrostatic  confinement  in  conjunction  with  the  full  line 
cusp  surface  magnetic  fields. 


2 


This  work  uses  a  standard  i+0  liter  multi -dipole  soup  pot  in  a 
full  line  cusp  configuration,^-  and  uses  magnetically  shielded 
electrodes  to  alter  the  internal  plasma  potential  structure  in  a 
way  that  improves  the  plasma  confinement. 

The  apparatus  used  for  this  work  is  described  in  Section  II . 
This  section  also  describes  briefly  some  of  the  diagnostic  devices 
used,  further  descriptions  being  provided  as  needed.  Section  II 
also  discusses  the  basic  operating  regime  of  the  device. 

Experimental  results  are  presented  in  Section  III.  Included 
are  variations  of  plasma  potential,  density,  and  temperature  with 
electrode  bias.  Also  presented  are  results  obtained  with  the  device 
configured  to  simulate  an  ion  source. 

A  simple  theory  is  presented  in  Section  IV  which  gives  good 
qualitative  agreement  with  the  results  of  Section  III. 

Section  V  deals  with  secondary  electron  emission  from  the 
chamber  wall  due  to  electron  bombardment.  The  appearance  of  the 
secondary  electrons  is  a  result  of  the  way  the  plasma  responds  to 
the  biased  electrodes  when  the  neutral  pressure  is  less  than 
10 Torr. 

Section  VT  is  devoted  to  seme  experiments  using  hydrogen  as 
the  neutral  gas.  Hydrogen  was  used  to  determine  the  effect  the 
electrodes  had.  on  the  ion  composition. 

Section  VTI  is  a  final  summary  of  the  work  presented. 


II .  APPARATUS 


The  chamber  used  for  this  work  was  a  40  liter  multi -dipole 
soup  pot,  with  the  permanent  magnets  in  a  full  line  cusp  arrange¬ 
ment.  The  magnets  produced  a  maximum  field  of  1  kilogauss  at  the 
poles . 

The  basic  device  was  modified  by  placing  positively  biased 
(with  respect  to  the  chamber  wall)  electrodes  straddling  each  row  of 

rr 

magnets  (Fig.  l) .  The  electrodes  were  approximately  0.7  cm  above 
the  wall.  At  this  distance  the  magnetic  field  was  of  order  5°0  Gauss, 
so  the  electrodes  were  magnetically  shielded  from  the  plasma.  It 
was  possible  to  change  the  gap  width  from  0.4  cm  to  1.3  cm.  This 
separation  is  on  the  order  of  the  ion-gyroradius  and  much  greater 
than  the  plasma  electron  gyroradius  (r^  s  1  cm  for  argon,  helium, 
and  hydrogen;  rg  s  5  X  10  ^  cm) . 

The  plasma  was  produced  by  electron  emission  from  negatively 
biased  filaments  located  in  the  field  free  interior.  The  energy  of 
these  ionizing  electrons  (hereafter  called  primaries)  was  equal  to 
the  difference  of  the  filament  bias  potential  cpp  and  the  plasma 
potential  cpp.  The  device  was  operated  with  low  electron  injection 
(discharge)  currents  (i^  s  1A)  which  produces  densities  on  the  order 
of  1010  cm”^,  and  low  fractional  ionization  (s  1%). 


The  primary  diagnostic  tools  used  were  collecting  and  emitt¬ 
ing  Langmuir  probes.  The  collecting  probes,  l/4"  diameter  Tantulum 
discs,  were  used  to  measure  plasma  electron  temperature  and  density, 
and  <Pp  in  the  field  free  regions.  The  emissive  probe  (3  mm  long, 
0.003  cm  diameter,  hot  tungsten  wire)  was  used  to  map  the  spatial 
plasma  potential  variation  in  the  device .  How  we  used  the  probes 
will  be  explained  more  in  Section  III  and  Section  V . 

This  work  also  required  the  use  of  mass  spectrometers  mounted 
inside  and  outside  the  chamber,  and  a  gridded  disc  for  launching 
pseudo-waves.  This  apparatus  will  be  discussed  more  fully  in 
Section  VT. 


Chamber  Operation 

The  primaries  are  confined  by  the  surface  magnetic  field. 
Previous  experiments  have  shown  that  the  leak  width  of  the  primaries 
is  only  the  order  of  their  gyroradius  (0.03  cm  for  the  maximum  1 B ! 

g 

of  1  kilogauss).  Plasma  electrons  and  ions  are  confined  by  a 
combination  of  magnetic  fields  and  self-consistent  electric  fields, 
which  for  large  electrode  gaps  (G  >  r^)  are  essentially  the  same  as 
those  in  multidipole  devices  without  electrodes. 

Consider  operation  with  the  electrodes  biased  at  such  that 
ecpp/T^  »  1.  For  large  gaps  the  plasma  potential  will  be  on  the 
order  of  Tg/e  above  ground,  which  is  the  anode.  As  the  gap  narrows 
to  less  than  the  self-consistent  ion  leak  width  the  ions  will  be 


electrostatically  confined  at  the  electrodes  and  the  ion  leak  width 


will  be  reduced  to  the  gap  width  or  slightly  narrower.  At  the  same 
time  the  electron  leak  rate  will  increase  because  of  the  presence  of 
the  positively  biased  electrodes.  Improved  ion  confinement  and 
reduced  electron  confinement  will  lead,  to  an  increase  in  the  plasma 
potential  until  the  electron  and  ion  loss  rates  balance.  As  the 
plasma  potential  increases,  electron  losses  to  the  wall  are  greatly 
reduced  (exponentially  as  e  ^p^e)  and  electron  losses  to  the  posi¬ 
tive  electrodes  are  increased.  Thus  the  electrodes  assume  the 
function  of  the  anode  in  the  plasma  and  the  plasma  potential  self- 
consistent  ly  adjusts  to  be  within  a  few  Tg/e  of  the  electrode 
potential. 

When  ***  <PE  >>  Te/e,  plasma  electrons  are  electrostatically 

confined  at  the  cusps,  as  mentioned  above,  and  plasma  electron  leaks 

are  then  dominated  by  diffusion  across  the  between  cusp  magnetic 

fields  to  the  electrodes.  Primaries  can  still  be  lost  at  the  cusps 

and  this  loss  flux  has  an  additional  consequence.  Cold  secondary 

electrons  with  energies  the  order  of  a  few  eV  are  produced  by  the 

primaries  at  the  walls  and  secondary  electron  currents  can  be  com- 

g 

parable  to  the  current  from  the  primaries .  In  ordinary  multidipole 
devices  these  secondary  electrons  are  indistinguishable  from  the 
plasma  electrons.  In  this  chamber  they  are  distinguishable  because 
the  secondaries  are  accelerated  through  e<p^/Te  »  1  before  reaching 
the  interior  of  the  device.  These  essentially  monoenergetic  electrons 
are  readily  confined  by  the  surface  magnetic  fields,  become  spatially 


isotropic  and  are  readily  identified  as  an  additional  electron 
species.  See  Section  V  for  more  details. 


The  energy  of  the  primaries  E  is  equal  to  e(cp  -  <p  )  because 

P  P  * 

the  plasma  serves  as  a  virtual  mode  for  the  filaments.  Since  the 


plasma  potential  cp^  generally  follows  (and  <PF  <  0),  cpp  was  varied 

to  keep  cp  -  <p _  =  constant.  The  ionizing  electron  energy  E  remains 
P  *  P 


essentially  constant  as  <Pg  is  varied.  Although  the  ionization  cross 
section  peaks  below  100  eV  (e.g.,  90  eV  for  argon);10  inelastic 


collisions,  multiple  ionizations  by  single  electrons,  and  other 
factors  combine  to  give  a  maximum  plasma  density  for  E^  greater  than 
100  eV,  and  this  optimum  E^  increase  with  increasing  neutral 
pressure.11  Above  this  energy  the  density  stays  close  to  the  maximum 


value.  For  this  work  the  primary  electron  energy  E  was  chosen 

P 

somewhat  greater  than  the  value  for  maximum  plasma  density  so  any 
residual  small  change  in  E^  does  not  have  the  effect  of  changing 
the  ionization  efficiency. 


III.  EXPERIMENTAL  RESULTS 


This  section  will  present  the  various  experimentally  deter¬ 
mined.  plasma  parameters,  plasma  potential  (cp  ) ,  electron  temperature 
(T  ),  and  density  (n)  and  their  variation  with  electrode  bias  (cp_), 
neutral  pressure  (pQ),  and  electrode  separation  (G) .  The  operation 
was  confined  to  two  separate  modes:  l)  raultidipole  plasma  chamber, 
2)  simulated  ion  source. 

The  diagnostics  used  were  emitting  and  collecting  Langmuir 
probes.  The  collecting  probe  was  typically  placed  near  the  center 
of  the  chamber.  A  voltage  was  swept  across  the  collecting  probe, 
and  an  I  vs.  V  Langmuir  probe  characteristic  was  generated  and 
recorded  on  an  x  -  y  plotter  (Fig.  2) .  The  method  of  getting 
parameters  from  these  characteristics  involved  removing  the  primary 
electron  portion  by  drawing  the  straight  line  through  that  part  of 
the  trace  (Fig.  3).  That  line  is  the  baseline,  the  zero  current 
line  for  plasma  electrons  (secondary  electrons  are  only  observed  at 
pressures  lower  than  those  used  in  this  part  of  the  experiment  and 
will  be  ignored  for  present),  and  the  natural  logarithm  of  the 
plasma  electron  current  is  plotted  against  the  applied  voltage 
(Fig.  k) .  This  semi -log  plot  will  ideally  have  two  parts.  First, 
a  linear  part  up  to  cp^,  and  a  second  part  that  is  not  collinear  with 


8 


the  first.  The  voltage  where  this  discontinuity  occurs  is  called 
the  plasma  potential.  The  voltage  difference  over  which  I  changes 
by  a  factor  e  is  the  electron  temperature  T  . 

There  is  a  third,  parameter  measured  on  the  characteristic  and 
that  is  the  plasma  saturation  current,  the  value  of  I  at  the  dis¬ 
continuity  mentioned  above.  This  current  can  be  written1"'" 


* 

I 

pe 


T 

n  A  (■— £ 


1/2 


2rr  m 


(1) 


A  is  the  probe  area,  and  Tg  is  measured  in  energy  units.  For  a 
l/4"  diameter  probe  the  density  can  be  written 


n  =  5.8  X  108  I^/JT  (2) 

-3  * 

where  n  is  in  cm  ,  I  is  in  mA,  and  T  is  in  eV.  Now  we  have 

pe  e 

measured,  or  calculated  the  plasma  parameters  from  the  character¬ 
istic. 


Multi -Pinole  Operation 

First  consider  the  variation  of  9^  with  applied  electrode 

voltage  9e  (Fig.  5)*  Note  that  cp^  generally  follows  q>E,  although 

initially  9p  >  9E;  then  after  e9p/TQ  >  2  the  opposite  is  true.  The 

difference  between  the  applied  voltage  and  plasma  potential  is  on  the 

order  of  T  . 

e 
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Tg  is  presented  (Fig.  6)  as  a  function  of  cpr.  Tg  is  shown  to 

be  a  fairly  constant  function  of  for  a  specific  parameter  set, 

the  error  bars  on  T  are  of  the  order  l/2  eV.  Later  it  will  be 
e  ' 

assumed  (Section  IV)  that  T  is  constant  and  these  data  show  that 

e 

this  assumption  is  warranted. 

The  density  also  increases  as  c?E  increases.  We  call  the  ratio 
of  the  density  to  its  value  with  zero  electrode  bias  the  enhancement 
ratio  (7).  Fig.  (7)  shows  the  variation  of  y  as  increases. 
Presented  also  in  this  figure  is  the  dependence  on  pQ  and  gap  width. 

Note  the  largest  ratio  of  y  is  obtained  when  G  =  12  mm  and 

-5  -L 

p  =  2  X  10  .  The  lowest  value  of  y  is  obtained  for  p  =  1  X  10 

o  o 

and  G  =  6  mm.  These  results  for  y  can  qualitatively  be  compared  to 
calculated  values  of  y.  This  is  done  in  Section  IV. 

Recall  this  experiment  was  initially  intended  to  prove  that 
the  electrodes  improve  the  ion  confinement,  but  ions  have  not  been 
discussed  yet.  We  infer  increased  ion  confinement  from  the  electron 
density  increase.  Physically,  what  is  impeding  the  ion  loss?  If 
the  ion  loss  area  is  reduced  then  the  electrodes  are  improving  the 
ion  confinement. 

Looking  for  the  reduction  in  loss  area  required  making  equi- 
plasma  potential  contours  around  the  cusp  area.  Since  this  area  was 
magnetic,  a  collecting  Langmuir  probe  couldn’t  be  used.  However, 

showed  that  an  emissive  Langmuir  probe  could  give  the 


Smith  et  al 
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desired  information.  We  used  the  inflection  point  method  they 
present. 

The  results  are  presented  in  Fig.  (8).  These  results  do 

show  the  formation  of  an  ion  trap.  In  this  case  the  ion  trap 

reaches  8  volts,  but  since  the  ions  axe  very  cold  a  trap  of  order  1 

volt  would  be  sufficient  in  these  devices. 

This  figure  also  shows  the  reduction  of  the  ion  loss  area. 

The  loss  area  is  defined  as  A_  =  W  L  where  L  is  the  length  of  the 

magnetic  cusp.  W^  is  an  empirical  loss  width  through  which  ions 

with  energy  e(«P  )  +  T  can  flow  out.  (<p  )  means  the  plasma 

p  c  1  p  c 

potential  in  the  center  of  the  device.  The  upper  figure  shows  this 

energy  is  E  =  0.3  +  T^  and  the  lower  figure  shows  =  h2  +  T^. 

The  corresponding  loss  widths  are  W  =  9  nun  and  W  -3  mm.  This 

L1  L2 

analysis  was  performed  for  many  <Pg,  and  a  linear  relationship  was 
found  between  n  and  (W^)  1  (Fig.  9) • 

Simulated  Ion  Source  Operation 

Operating  the  chamber  as  a  simulated  ion  source  was  done  by 

covering  one  end  of  the  cylinder  by  a  grid  (or  a  plate)  and  removing 

the  magnets  and  electrodes.  We  again  measured  5p^,  n,  and  Tg  as 

functions  of  pQ.  Gap  width  was  set  at  12  mm. 

Representative  data  are  shown  in  Fig.  (10).  The  points  to 

notice  are:  l)  V  is  slightly  less  than  «p_  at  high  electrode  bias 
P 

and  insensitive  to  whether  the  grid  is  floated  or  grounded,  2)  the 


density  does  saturate,  but  the  saturation  occurs  at  higher  cpr  than 
when  the  grid  is  not  present,  and  3)  the  factor  by  which  the  density 
increases  is  greater  for  the  grounded  grid  case.  A  simple  calcula¬ 
tion  explaining  the  variation  of  density  for  a  grounded  and  floating 
grid  is  presented  in  the  next  section. 

Tg,  for  the  case  shown  in  Fig.  (10),  is  graphed  in  Fig.  (ll). 
Again,  Tg  is  a  fairly  constant  function  of  q>v. 

Thus,  overall  the  grid  has  very  little  impact  on  the  varia¬ 
tion  of  the  measured  parameters  of  the  plasma  with  electrode  bias. 

Of  course,  the  density  is  reduced  because  the  grid  presents  a  large 
loss  area  for  ions  and  electrons,  there  being  no  magnetic  shielding 
of  the  grid.  However,  the  main  point  is  that  whether  the  chamber  is 
operated  as  an  ion-source  or  multidipole  soup  pot,  the  results  are 
qualitatively  the  same  (compare  Fig.  7  and  Fig.  10) . 


IV.  THEORY 


Numerical  Estimate 

This  section  will  present  a  simple  model  that  predicts  a 
density  increase  in  qualitative  agreement  with  the  experiment.  A 
basic  assumption  used  at  this  time  is  that  the  plasma  potential  is 
increased  as  the  electrode  bias  is  increased,  and  about  equal.  In 
the  second  part  of  section  IV  this  assumption  will  be  proved. 

The  quantities  to  consider  are  the  electron  loss  rate  along 
the  cusp  at  low  electrode  bias,  and  the  electron  loss  rate  across 
the  magnetic  field  at  high  electrode  bias.  These  quantities  are 
considered  because  they  are  the  dominant  electron  loss  mechanisms. 

The  electron  loss  rate  out  the  cusp  can  be  expressed  as: 


n  v 
e  e 


eLc 


w  L  e 
c  c 


-  e  <p  /T 
p'  e 


(3) 


The  ion  loss  rate  can  be  written: 


vllc  -  f  (AIC  *  Afll' 


(>1) 


j+/e  is  a  Bohm  loss  flux  which  equals  ncg/2.  is  the  ion  loss 

area  at  the  cusps  and  equals  p  .  P  is  an  effective  ion  leak  which 
has  been  shown ^  for  similar  devices  to  equal  Vrg  r^,  Lc  is  the 
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total  cusp  length.  the  area  of  the  filaments,  is  only  0.1$ 

of  so  we  will  ignore  it.  ng  is  the  total  primary  and  plasma 
electron  density,  but  the  primary  density  is  small  so  we  approximate 
ng  by  n,  the  interior  plasma  density.  vq  is  the  electron  thermal 
velocity  (v  =  V(Tg/M  ) ,  and  cg  is  the  ion-acoustic  velocity 
(c  =  J'l JV[. ) .  The  exponential  takes  into  account  electron  electro- 
static  confinement  along  the  cusp.  The  only  quantity  undefined  is 
w  which  is  an  effective  electron  loss  width,  and  the  unknown  of  the 
present  calculation. 

Equating  the  electron  and  ion  loss  rates,  because  charge 
neutrality  must  be  maintained,  allows  expression  of  w^  in  terms  of 
known  parameters.  Therefore, 


t  v  w  L  e 
4  e  c  c 


e  4  /T  nc 

«rV  a. 


*  e  «  -J-  Vr“r7  L 
2  e  i  c 


c  e  «P  /T 

o  s  r—  p  e 

w  =  o  — vr  r.  e  v 
c  v  e  i 
e 


1/4 


m  T.  X/'T  e  cp  /T 

TN'  C 

e 


.8,  (-fA  •  V  e 


e  m.  T 
i  e 


(5) 


(6) 


Assuming  an  argon  plasma  with  Tg/T^  =  10  gives 

e  tp  /T 

_  ^  tv  e 

w=  0.27  r  e  F 
c  e 


(7) 


lb 


Of  course,  w is  determined  by  physics  other  than  that  involving  the 
interior  plasma  potential  cp^.  In  fact,  if  is  determined,  then 
this  equation  determines  the  plasma  potential  such  that  the  electron 
and  ion  loss  rates  are  balanced.  For  example,  using  wc  =  rg,  then 
gives  e  (P  *  1.29  Tg,  a  result  which  is  very  close  to  that  experi¬ 
mentally  observed  in  the  present  work.  Thus,  w^  cannot  be  greatly 
different  from  rg  or  the  predicted  plasma  potential  will  not  agree 
with  the  data.  Furthermore,  in  a  discussion  of  the  properties  of 
the  systems  as  the  electrode  bias  is  changed,  the  electron  loss 
along  the  cusp  will  be  completely  dominated  by  the  exponential  terms 
in  Eq.  (l)  so  the  exact  value  used  for  w^  is  unimportant.  Therefore, 
we  take  wc  =  rg  in  the  following  and  write 


v 

e 


r 

e 


L 

c 


4  V 


(3) 


Now  consider  the  case  e  V  /T  »  1.  This,  of  course,  results 

P  ® 

from  a  high  electrode  bias.  In  this  regime  the  electrons  are 
electrostatically  trapped  along  the  cusp,  so  the  only  loss  process 
for  electrons  is  diffusion  across  the  magnetic  field  lines. 

Dolan  defines  a  characteristic  time  for  electron  loss  by 
spatial  diffusion  to  be:^ 
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T  =  (total  number  of  electrons) 

d  _  (number  lost  per  second  by  diffusion) 

=  -  1  g.  n _  (9) 

j  d4  •  (-  D  V  n) 


■where  is  a  volume  element  and  di£  is  a  surface  element  of  the 

plasma.  The  surface  element  can  be  written  da  =  Ldl  where  L  is  the 
height  of  the  plasma  cylinder  and  dl  is  along  the  magnetic  field. 
.'Jew 


dl  n  *  r.  L  tt  r  c 
P 


(10) 


with  r 

P 


the  radius  of  a  cross  section  of  plasma. 

Tk 

The  diffusion  coefficient  is  given  by: 


D 


m  T 


e^2 


v 

c 


(11) 


where  is  the  effective  momentum-transfer  collision  frequency  for 
electrons  with  ions,  neutrals,  and  E  field  fluctuations  in  the  limit 
vc  «  (u^^.  The  density  gradient  may  be  approximated: 


-  Vn  »  n/^  (12) 

where  is  the  distance  between  flux  surfaces  and  ^-(Fig.  12). 
is  the  flux  surface  that  just  touches  the  inner  edge  of  the 
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electrode,  and  ♦.  is  a  flux  surface  that  is  a  distance  r  frcm  the 
’  0  e 

center  of  the  magnetic  cusp. 

Plugging  these  factors  into  Eq.  (9)  gives 


T 

d 


2  2 

it  rp  e  /me 

P  cU  T  v  /r  32 
J  e  c/  1 


(13) 


Let  ?0  and  BQ  be  the  values  of  5^  and  3  at  any  arbitrary  reference 
point  (Fig.  12).  By  magnetic  flux  conservation: 


-Vo 


(lM 


v  and  T  are  constant  in  the  main  nlasma  volume  because  the 
c  e 

central  plasma  is  field-free  and  uniform.  We  will  assume  that 

and  T  are  constant  all  along  By  definition  dl#  =  di/2rr  r  and 

el  p 

B'  =  B/B  .  Then: 


e2  B  2  l  r 

T _ 2  0  p  _ 

d  2m  T  v  X 


J&Ie— 

2 

I  V  X 
o  c 


(15) 


X  is  written:1^ 


X 


«  1 


(16) 
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N  is  the  number  of  magnetic  cusps,  and  PQ  is  the  electron  gyro-radius 

at  the  reference  point  where  5^  and  B  are  evaluated. 

o  o 

The  cross  field  leak  rate  for  electrons  is  defined 


where  V  is  the  plasma  volume.  Now  substituting  gives: 


p  2  v 

',«cr =  nV 


-  s  e(«Pp  -  <PE)/Te 


Now  V/r  =  A  the  surface  area  of  the  plasma: 


A  «  SL  , 
c  * 


with  S  the  distance  between  cusps  (**  8  cm)  and  Lc  the  magnetic  cusp 
length.  Now  5Q  is  the  order  of  half  the  electrode  gap  G.  The 
electron  gyro-radius  at  the  boundary  ^  can  be  written  PQ  =  vg/ft 
where  fi  is  electron  cyclotron  frequency.  The  exponential  is  again 
an  electrostatic  confinement  factor  with 


CPp>5PE 


'PE>‘Pp 
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N  is  the  number  of  magnetic  cusps,  and  pQ  is  the  electron  gyro-radius 
at  the  reference  point  where  5q  and  Bq  axe  evaluated. 

The  cross  field  leak  rate  for  electrons  is  defined 


v  =  & 
eLT  r 

a 


(17) 


where  V  is  the  plasma  volume.  Now  substituting  t A  gives: 


v  .  „v  .  '  8  ■  %)/Te 

?orp 


(18) 


Now  V/r 


P 


=  A  the  surface  area  of  the  plasma: 


A  =  SL 

c 


(19) 


with  S  the  distance  between  cusps  (*»  8  cm)  and  the  magnetic  cusp 
length.  Now  is  the  order  of  half  the  electrode  gap  G.  The 
electron  gyro-radius  at  the  boundary  can  be  written  PQ  =  vg/ft 
where  ft  is  electron  cyclotron  frequency.  The  exponential  is  again 
an  electrostatic  confinement  factor  with 


g  = 


fp  >’’E 


cp_  >  m 
^E  p 
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The  cross  field  leak  rate  is: 


e.  .  _  ”  '  'p 


veLT  -  n  ve  re  Lc  ‘  (TT  '  T)e 


-  g  e(cp  -  <?J/T 


(20) 


_  n 


r  v  w  L  e 
4  e  T  c 


g  e(cpp  -  «PE)/Te 


(21) 


Where  in  the  last  equation  we  define  an  effective  cross -field 
leak  width,  wT,  given  by 


,  c  2S 

WT  =  Ure  T  T  * 


(22) 


can  be  written,  ignoring  E  field  fluctuations,  as 


v  =  v  .  +  \> 
e  ei  en 


(23) 


where  v  .  is  the  electron-ion  collision  frequency  and  ^  is  the 
ei  ^  en 

electron -neutral  collision  frequency,  v  ,  and  v  are  estimated 

J  ei  en 


as : 


.15 


\>ei  =  4  X  10  '  n/T 


-5  5/2 


v  =  5.5  x  108  p  T  5/ 2 
en  *o  e 


(24) 

(25) 
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The  density  n  is  measured  in  cm  ^ ,  pQ  in  Torr,  T^  in  eV,  and  \> 

-1 

in  sec 

Q  is  the  cyclotron  frequency  evaluated  midway  between  the 
cusps.  Measuring  the  field  strength  along  it  was  found  that  | B | 
at  this  point  is  approximately  20  Gauss.  Taking  |b|  =  0  in  the 
interior,  we  see  an  average  |b|  for  electron  diffusion  is  10  Gauss. 

g 

This  gives  Q  =  1.8  X  10  rad/sec. 

Now  taking  the  ratio  of  wc  and  w^,  assuming  Te  ~  3  eV  gives 


w  r 

_c _ e_ 


w_ 


4r  v  __ 
e  c  2S 

n  G 


G  fl 
8S  v 


=  (2.8  X  10°  G) 


x  (7.7 


-6 

X  10  n 


+  2.86  X  10' 


*o> 


(26) 


This  expression  with  WC/WT  normalized  to  G  is  graphed  in  Fig.  (l3). 

-4  9  -3  10  -3 

This  graph  shows  that  for  pQ  ~  10  and  10^  cm  *  n  *  10  cm 

then  (w  /w„)/G  ~  8.  Since  0.6  cm  £  G  s  1.2  cm  we  see  that 

4.8  *  w  /v-  £  9.6. 
c'  T 

In  general,  plasma  density  is  determined  by  a  balance  of  plasma 
production  and  losses.  The  creation  rate  of  plasma  electrons  can  be 
written  as  81^  where  1^  is  the  filament  discharge  current  and  6  de¬ 
pends  on  the  cross  section  for  ionizing  and  nonionizing  collisions, 
and  on  the  primary  electron  confinement  time.  During  an  individual 
experiment  81^  was  constant,  independent  of  cp£,  therefore 
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eId  = 


eLc 


+  v 


eLT 


(27^ 


When  cp_  =  0  and  e  fp  /T  ~  1  means  31,  ~  s>  _  ,  but  when  9-  ~  <P  »  T  /e 
E  p'  e  d  eLc7  E  p  e' 

(cpg  ~  cp_^  will  be  proved  later)  then  31^  ~  •'-as't  sentence 

refers  only  to  the  dominant  loss  mechanism  at  a  given  limit  for  cp^. 
Equating  the  two  loss  rates  shows 


L 

c 


-  e  cp  (0)/T 
p  e 


T  ve  WT 


L 

c 


(28) 


The  density  ratio  is 


w  -  e  cp  (0)/T 

_£.e  P  e 

wm 


w 

0 .27  f 
WT 


(29) 


Using  the  limits  of  wc/w^,  means 

1.25  s  — s  2.5  .  (30) 

n 

o 

These  values  are  in  excellent  agreement  with  the  ratios  observed  in 
Fig.  (7). 

At  sufficiently  high  neutral  pressure  or  low  plasma  density 

where  electron-neutral  collisions  dominate  (v  _  h  __)  a  G/v  aG/p  . 

eLc'  eLi  en  o 

However,  if  we  consider  the  case  of  v  »  v  then 
’  ei  en 


2] 


a  a  - 

veLT  vei  n 


(31) 


Here  the  effective  cross  field  leak  rate  is  proportional  to  the 

plasma  density  so  the  steady  state  density  ratio  is  no  longer 

proportional  to  the  ratio  given  by  Eq.  (29).  The  neutral  pressure 

and  density  at  which  the  transition  from  v  »  v  to  v  .  »  v 

en  ei  ei  en 

depends  on  the  particular  plasma  species  and  electron  temperature. 
The  dependence  of  the  enhancement  ratio  on  the  gap  width  G  and  the 
discharge  current  1^  can  he  seen  as  follows.  At  high  electrode  bias 
the  balance  of  electron  production  and  loss  (when  v  »  v  )  can 
be  expressed: 


Hd  “  "H  vr  r 


-  L  T- 

e  e  c  G 


(32) 


so  the  density  at  high  electrode  bias  goes  like  n^  a  Vl^G  .  At  low 
electrode  bias,  when  losses  along  the  cusp  dominate,  the  density  is 
proportional  to  1^  so  the  ratio  of  density  at  high  electrode  bias  to 


zero  bias  is: 


a.  vCG 

n.  1A 
o  d  Td 


(33) 


in  agreement  with  Fig.  (7)' 


:  vlTi 
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Variation  of  n  with  q>E  In  a 
Multi -Pinole  Field 

The  dependence  of  the  density  on  electrode  bias  can  be  deter¬ 
mined  from: 


v 

e 


L  [w 
c  c 


-  e  cp  /T 

ty  e 

e  v  +  wT  e 


(34a) 


However,  Eq.  (34a)  requires  a  relationship  between  tp  and  cp^.  To 
find  this  relationship  requires  equating  electron  production  and  loss 


£  v  L  [w  e 
4  e  c  c 


e  cp  /T 
p7  e 


*  3  e(l>  -  <Pr)/T 


+  e 


!] 


c 


s 


Vr  r.  L 
e  x 


c 


(34b) 


where  R  is  a  retardation  factor  that  is  unknown,  and  v  »  m  is 

en  ei 

assumed.  This  equation  can  be  reduced  to 


-  e 


V*.  ,1  ■  8  e(Tp  -  %>/Te 

^  +  O  €  =  Ue 


27  R 


(35) 


where  w  /w_  ~  8  has  been  used, 
o'  T 

Consider  first  <p_  =  0,  but  9  >0  then  Eq.  (35)  can  be  written 
“  P 


1 

8 


<P  /T 

V  e 


0.27 


(36) 


25 

where  R  =  1  because  the  ions  are  lost  by  an  unin^eded  Bohm  flux. 
Therefore 

e  cpp/Te  -  1.4  ■  e  cpp(G)/Te  . 

Consider  now  <p„  >  0,  but  cp  =  cp  (o).  When  tp  =  cp  (o)  Eq.  (57) 
b  P  P  bp 

becomes 


|  (0.27)  +  |  =  0.27  R  =  0.5678 


(38) 


or  implies  R  >  1.  However,  R  ^  1  so  cpp  must  increase  as  some 
function  of  cp^,  to  keep  R^l.  Since  this  function  is  unknown  R 
cannot  be  calculated,  but  limits  of  R  can  be  obtained. 

Assume  initially  that  R  =  1  and  <pp  >  <PE*  Therefore  Eq.  (39) 
is,  solving  for  <pp: 


+  i 


e  cp  /T  =  2n(- 
p'  e 


0.27 


(40) 


This  equation  satisfies  cpp  >  cp^,  only  if  e  5PE/Tg  ^  2.  Therefore,  R 

must  be  less  than  1  when  e  <PE/Te  >  2.  Let  us  now  look  for  R,  such 

that  <10,  =  cp  and  define  this  R  =  R  .  Therefore,  Eq.  (55)  becomes 
E  p  max 
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-  e  <Pe/T 


e 


0.27 


R 

max 


(41) 


*/ 6.4 
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"  e  {Pe/t^ 

or  R  =  3. 7(l/8  +  e  ).  For  this  experiment  0  *  e  tp„/T 

ni &yz  ^  6 

s  20  (T  ~  3  eV),  therefore,  0.46  s  R  £  1.0. 

Consider  sp^  <  q>E,  now  the  only  relationship  between  <p^  and  cp^, 
is  through  R.  This  regime  is  the  case  considered  in  Eq.  (40) 
with  <p^  replacing  t?E.  That  is,  for  R  as  large  as  is  physically 
possible,  cp^  =  cp^ .  We  conclude  that  even  if  ion  losses  are  retarded 
only  to  the  extent  that  they  can  physically  balance  electron  losses, 
we  must  have  cp^  ~  q>  . 

The  smallest  R  that  is  physically  reasonable  occurs  when  the 
cusps  are  plugged  by  the  electrode  potential  then 


(42) 


where 


*E  >CPp 


h  =  < 

0 


Eq.  (35)  can  be  written 


-  e  cp  /T  .  ■ 

e  *  *  *  J  e 


8  e(fp  -  »E)/Te 


=  0.27  e 


-  he(«E  -  •  )ft1 


(43) 


Note  if  g  =  1  then  h  =  0  and  vice  versa.  There  cannot  be  plugging 
for  electrons  and  ions  at  the  same  time.  Eq.  (43)  can  be  solved 


numericaUy  for  cp^  as  a  function  of  cp^.  The  solution  set  of  Eq.  (43), 
for  Tg/^  =  10  and  <  c?E,  is  that  ete^  -  ep  )  <  Tg  and  e(cp£-  cp^)  is 
largest  when  e  qJg/Tg  ~  20  [i.e.,  e(«pE  -  cp^)  =  0.08  Tg] .  So  the  plasma 
potential  is  slightly  less  than  the  electrode  potential.  We  conclude 
from  this  that  at  elevated  values  of  cpE,  ~  cp£  and  we  will  assume 
equality  in  the  following. 

Let  us  now  consider  the  density  dependence  on  tp_.  Using 
w^/wc  =  1/8,  and  <PE  =  <p^  »  0,  the  balance  of  electron  production 
and  loss  gives 


pi ,  =  r  v  L  w  [e 
d  4  e  c  c 


e  %nc 


(44) 


When 


=  0  then 


v  L 
e  c 


w(l 
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-  e 


cp  (0)/T 
p  e 


(45) 


The  density  ratio  is 


"  6  'Pr)(0)/Te 

n  2  e  p  e 
n(o)  ~  8  -e  %/T  * 

e  E  e  +  1/8 


This  equation  is  graphed  in  Fig.  (13). 

Now  assume  »  vgn  then  Eq.  (44)  becomes 
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where  cp^. 


«  cp 

P 


»  0. 


Then: 


n  2L  ® 
n(0)  ”  8 


-  e  cp  (0)/T 
P  e 

'  e  VTe  _jl_ 

8n(o) 


(U8) 


This  expression  is  also  graphed  in  Fig.  ( 13) . 

The  above  considerations  show  that  a  theory  assuming  the  main 
loss  mechanism  for  electrons  at  large  cp„  is  cross  field  diffusion 
gives  very  good  qualitative  agreement  with  the  experimental  results 
of  Section  III. 


Variation  of  n  with  cpr  in  the  Presence 
Pf ..an.Fxtrftp.tjon  r,Fj d 

Now  consider  a  problem  similar  to  the  one  just  discussed,  but 
where  one  end  of  e-pot  is  replaced  by  a  grid  or  metal  plate.  The 
calculations  discussed  here  will  show  good  qualitative  agreement 
with  experimental  results  presented  earlier  about  simulated  ion- 
source  operation. 

First,  consider  the  potential  of  a  floating  grid.  The  ion 
current  to  such  a  grid  with  area  A  and  at  a  potential  cp  can  be 

e  g 

estimated  to  be  I  =enc  A/2  provided  cp  >  cp  .  T  will  only  have 
g  ®  g  P  &  g 

variation  with  n,  as  we  showed  earlier  T  is  constant,  and  A  is 

e  g 

constant.  When  the  grid  potential  V  =  the  plasma  floating  po- 

£  * 

tential,  the  plasma  must  supply  electron  current  to  balance  the  ion 
current.  At  high  neutral  pressure  where  the  density  of  the  primary 
electrons  is  much  less  than  the  plasma  density  this  balance  gives 
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n 


-e(»p  -  <PfVTe  _  a 


A  J  e 


or  solving  for  <p  -  <P„  gives 
p  f 


V  -  <p V  =  «  5  T  /e 

P  re-  —  e' 

*  m  2rr 

e 


for  an  argon  plasma.  However,  at  low  neutral  pressure  the  primary 

g 

electrons,  with  energy  E  ,  contribute  a  current  1^  given  by: 


I  =  e  n  A  v  cos  9 
P  P  g 


where 


v  cos  6  =  J  d$  fO$)v  cos  9 


The  distribution  function  fC$)  is  characterized  as  a  shell  in 

,  ..  16-17 

velocity  space 


ftt)  = 


where  v^  is  primary  electron  velocity.  Therefore 


•MV****.  1'  •  -  .  /»rj 


where  eV  =  e(cp  -  cp  ) .  Then  the  expression  for  I  is 
PS  p 


v  E  -  e  (q>  -  cp  ) 

I  =  e  n  -5*  A  (-2 - n2 - 

P  P  U  S  E 


(55: 


e  n  _  E  -  e(cp  -  cp  ) 

-S*  VEE  (T  eT  *> 
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The  current  balance  gives 


(56) 
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E  -  e(cp  -  cp  ) 
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AJTTT  ( J - —2 - S_)  +  e  n  aJ~T! 


-  e(cp  -  cp  )/T 
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S  _e_ 

2na 


P  S  '  e 


§Vt  a 
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(57) 


Since  the  right  hand  side  of  Eq.  (49)  is  the  same  as  Eq.  (57) 

requires  that  cp  -  cp  must  increase  (reducing  the  plasma  electron 
P  S 

current).  This  trend  in  cp  -  cp  is  shown  in  Fig.  (l4). 

P  g 

A  similar  calculation  as  was  used  earlier  this  section  can 
employed  to  calculate  the  density  variation  with  cp^  when  the  grid 
is  in  place.  With  a  floating  grid  the  balance  of  electron  produc¬ 


tion  and  loss  can  be  written 


! 


v  L  ■  e  <J  /T  A 

0Id  -fT^V  P  6  +  V  +  Vsf 


-  e<p  /T 

nf  a  [wc  e  P  6  +  wT  +  2  cs  Ag/v0  L.]- 


This  expression  ignores  primary  electrons  so  is  only  valid  at  high 


pressure.  This  equation  also  shows  that  density  is  reduced  by  the 


presence  of  the  grid,  but  the  variation  with  cp£  is  not  qualitatively 


changed. 


If  the  grid  is  grounded,  in  the  same  limit  as  the  floating 


grid  case,  the  density  varies  as: 


■  e  T  /T 


n  a  (w  e  p  e  +  w„  +  'A  e 


gn  c 


4  A  -  e  cp  /T 

r  +  - s.  e  P'  *) 

T  L  ’ 

c 


When  cp  <0  electron  current  will  be  less  than  when  cp  =  0.  Then 
8  g 


the  last  term  in  the  expression  for  n^,  will  be  less  than  for  the 


last  term  in  Eq.  (60) .  Then  the  ratio  of  (n„/n  )  <  1.  If  cp  >  0 

x  gn  g 


the  opposite  holds  and  n^/n^  >  1.  These  arguments  are  in  quali¬ 


tative  agreement  with  the  data  shown  earlier  in  Fig.  (10). 
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V.  SECONDARY  ELECTRONS 

One  consequence  of  the  elevated,  interior  plasma  potential  is 
the  observation  of  secondary  electrons  by  primary  electron  bombard¬ 
ment  of  the  chamber  walls.  These  secondary  electrons  have  a 
characteristic  signature  in  the  Langmuir  probe  trace  that  is  clearly 
distinct  from  the  plasma  and  primary  electron  contributions. 

The  physical  picture  is  a  cold  electron  produced  at  the  wall 

being  accelerated  through  the  plasma-wall  sheath  at  the  gap  and 

forming  a  shell  in  velocity  space  at  v  =  «/2  e  cp  /M  .  In  ordinary 

P  e 

multi -dipole  devices  this  velocity  is  of  order  ./Te/M^  and  the 
secondaries  merely  spread  out  the  Maxwellian  plasma.  The  signature 
left  by  secondary  electrons,  when  projected  onto  a  planar  probe,  is 
a  straight  line  starting  at  tp  =  0  (Fig.  15 ).  This  straight  line  is 
most  easily  seen  when  Pq  <  1  X  10  .  At  the  lower  pressure,  the 

density  of  primaries,  and  therefore  secondaries,  is  a  large  fraction 
of  the  total  electron  density  and  therefore  easily  observed.  Also, 
at  these  low  pressures  the  secondaries  do  not  readily  thermalize 
with  the  plasma  electrons. 

The  number  of  secondaries  produced  per  incident  primary  is 
called  the  yield  (6).  The  6  wanted  is  for  production  at  the  wall 
and  should  not  include  production  at  the  probe. 


The  probe  current  can  be  written 


I  * 

v 


vf (v) [l  -  S(v  ) ]dv 
w 


where  8 (v  )  is  the  ratio  of  secondary  current  to  primary  current. 

W 

The  velocity  of  electrons  at  the  probe,  v^,  goes  to  zero  when 

2  2 

v  =  v  .  .  Energy  conservation  gives  v  +  2  e  tp  /m  =  v  so 
min  pr  w 

v  .  =  */  -  2  e  cp  Tin"  where  is  the  urobe  bias  and  e/m  is  electron 

min  pr'  pr  *  ' 

charge  to  mass  ratio.  f(v)  is  the  velocity  distribution  for  a  shell 

IS 

of  monoenergetic  electrons  projected  onto  a  one  dimensional  probe 


0  s  v  s  v 


v  <  v 
P 


where  v  is  the  primary  electron  velocity.  6(v  )  is  known  to  vary 
p  w 

with  bombarding  energy  as  shown  in  Fig.  l6.  For  simplicity  we  assume 


6  (v  )  »  ( — ) 
w  v 

u 


where  v  is  defined  by  =  l/2  My^  ,  and  E^  is  the  secondary 
electron  energy  such  that  8=1  (Fig.  1 6).  Thus,  the  probe  current  is 


V  c:  / 

p  v  +  2  e  cp  /m 

I  «  J*  v(l  -  - r-pr— )  dv 


or 


1  *  i-  f (E„  * e  %)(Eu  - e  v’  ‘  K2 '  '2  ,-2)1  (65) 


pr 


■where  E^  =  (l/2)  mv  .  The  primary  electron  current  to  the  probe 

when  the  probe  is  biased  at  <P  ,  is  the  saturation  current 

I  oe  (e  +  e  cp  )/m. 

P  P  P 
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Hachenberg  and  Brauer  list  E^  89  200  V  for  tantalum,  which 

the  probe  was  made  from.  The  assumed  form  for  8(v  )  is  a  slight 

w 

underestimation,  but  the  maximum  error  is  about  25%  at  E  =70  eV. 

P 

A  representative  probe  trace  which  exhibits  all  three  electron  species 
and  effects  from  the  probe  is  shown  in  Fig.  (l6).  Note  the  satura¬ 
tion  currents  of  the  three  components  are  comparable  and  that  the 
reduction  in  the  apparent  primary  electron  current  is  significant. 

The  fit  of  the  three  components  is  not  sensitive  to  the  yield 
coefficient  used  for  secondary  electrons  emitted  at  the  probe. 

Now  that  the  saturation  currents  have  been  corrected  they  can 
be  used  to  calculate  the  yield  coefficient  at  the  chamber  walls. 

Secondaries  are  created  at  a  rate  given  by  6n  v  A  where  n  is  the 

P  P  P  P 

density  of  primaries,  v^  is  primary  electron  velocity,  and  is 

nrimary  loss  area.  Secondaries  are  lost  either  to  the  walls  (A  )  or 
*  s 

by  inelastic  processes  identical  to  those  involving  primaries.  If 
the  secondaries  cause  an  ionization  they  lose  enough  energy  to  be 
trapped  and  then  lost  by  diffusion  (recall  Section  IV).  Balancing 
production  and  loss  gives 
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5n  v  A  =  n  v  (A  +  Vo  n  ) 
ppp  sss  io 


where  V  is  the  plasma  volume,  nQ  is  the  neutral  density  and  is 
the  ionization  cross  section.  The  probe  saturation  currents  for 
primaries  and  secondaries,  I  ,  are  proportional  to  nv  so  Eq.  (66) 
gives 


(67) 


The  term  A  =  r  L  where  r  is  the  primary  electron  gyroradius , 
ppc  p 

and  L  is  total  length  of  the  magnetic  cusps.  Now  r  =  2.V2E  /B  . 
c  P  P 

Recall  from  Section  II  that  at  the  wall  B  ~  1  kilogauss  so 

r  OA  X  10‘3  */E~.  A  typical  energy  for  primaries  is  E  ~  100  eV 
P  P  P 

2 

or  A  =  27  cm  . 

P 

The  numerator  of  the  bracketed  term  is  not  as  easily  calcu¬ 
lated.  It  is  possible  for  A  to  be  as  small  as  r  L  or  as  large 

s  sc 

as  GL  where  r  is  the  secondary  electron  gyroradius  and  G  is  the 
c  s 

2  2 

electrode  gap,  so  6  cm  s  Ag  £  800  cm  .  The  ionization  term,  for 

an  argon  plasma  at  neutral  pressure  pQ,  can  be  evaluated  using 

V  =  2.2  X  104  cm5,  <7,  =  2.4  X  lO-16  cm2,  and  n  =  3  X  10l6  p  (Torr) 
'  i  o  o 

or 

V<7,  n  »  1.6  X  105  p  (Torr) 
i  o  o 


(68) 


-6  o 

At  p  =  5  X  10  Torr,  this  term  is  .8  cnT.  Therefore  the  ioniza- 
o 

2  2 
tion  is  negligible  if  A  =  800  cm  and  non-negligible  if  A  =  8  cm  . 

s  s 

The  extreme  estimates  of  the  bracketed  term  gives  0.33  and  30.0. 

Although  E^  varied  for  data  in  this  experiment,  E^  is  determined  by 

the  details  of  secondary  emission  and  essentially  constant  so  the 

bracketed  term  can  be  written  aA/E-  . 

P 

Experimental  values  of  I  /I  are  shown  in  Fig.  (17).  The 

s  p 

*  * 

values  of  I  /I  are  comparable  to  1.0  so  the  minimum  value  of  6  for 

s'  p 

E  =  100  eV  is  0.33*  A  graph  (Fig.  17)  of  6  =  (i  /i  )  IOA/eT"  shows 

p  s'  p  '  p 

the  general  initial  trend  for  6 ,  where  10  is  a  representative  value 

taken  for  the  numerator  of  the  bracketed  term  in  Eq.  (67). 

The  above  argument  has  shown  that  primaries  can  produce  a 

significant  number  of  secondaries.  However,  the  main  point  was  to 

consider  the  secondary  electron  contribution  to  the  plasma.  The 

production  rate  of  secondaries  is  6n  v  A  while  the  production  of 

P  P  P 

plasma  electrons  is  n^  v^  nQ  V,  ignoring  primary  and  secondary 
electrons  degraded  in  energy  by  collisions.  The  ratio  of  secondary 

r 

production  to  plasma  production  is  S  A^/V  <7^  nQ  =  3^  6  when  A  =  27  cm£ 
and  pQ  =  5  X  10  ^  Torr.  Using  the  minimum  of  0.33  for  6  means  that 
6  A  nQ  ^  10  so  more  electrons  are  being  produced  by  secondary 

emission  than  by  ionization. 

It  has  been  shown  that  plasma  electrons  are  electrostatically 
confined  at  the  cusps  and  losses  occur  by  diffusion  between  the 
cusps.  The  secondaries  are  not  electrostatically  confined  'unless 


they  lose  energy.  The  ratio  of  density  of  secondary  electrons  to 
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plasma  electrons  can  be  quite  high  and  their  probe  saturation 
currents  can  be  comparable.  Since  in  conventional  multidipole 
devices  e  **  1,  the  confinement  of  plasma  and  secondary 

electrons  is  essentially  the  same  and  this  means  neglecting  the 
secondary  electrons  could  introduce  large  errors  in  calculations 
involving  plasma  production  and  loss  rates. 

Secondary  electrons  have  not  been  considered  in  previous 
multidipole  experiments  because  they  are  indistinguishable  from 
plasma  electrons  -when  wall  sheaths  are  small.  In  addition,  at  high 
neutral  pressures  primary  electron  densities  are  relatively  low. 

The  presence  of  secondary  electrons  is  not  negligible  for  low  neutral 
pressure  regimes.  Calculations  of  particle  and  energy  balance  in 
multidipole  devices  that  do  not  take  secondaries  into  account  can 
result  in  significant  errors.  Secondary  electrons  are  also  impor¬ 
tant  in  interpreting  Langmuir  probe  data.  It  is  apparent  that  the 
secondary  electrons  emitted  from  the  probe  can  lead  to  substantial 
misjudgement  of  the  primary  electron  saturation  current  and  conse¬ 
quently  the  plasma  electron  temperature  (Fig.  L). 


VI.  RELATIVE  ION  DENSITIES 


This  section  deals  with  experiments  where  the  relative  density 
of  the  various  ions  species  are  observed.  Directly  measuring  the 
relative  ion  densities  would  give  clear  evidence  concerning  effec¬ 
tiveness  of  electrodes  in  confining  the  ions .  Showing  that  the 
atomic  ion  density  increases  with  electrode  bias  relative  to  the 

molecular  ion  density  could  have  a  large  effect  on  whether  ion 
20 

sources  for  neutral  beam  injectors  are  modified  with  electrostatic 
plugging.  Unfortunately,  this  measurement  could  not  be  performed. 

The  relative  density  measurement  was  attempted  by  three  methods: 
l)  mass  spectrometer  mounted  outside  the  chamber,  2)  particle  bursts 
called  pseudo-waves,  and  3)  mass  spectrometer  mounted  inside  the 
plasma  volume. 


Theory 

The  idea  behind  a  mass  spectrometer  is  fairly  well  known,  but 
to  refresh  the  reader's  memory  the  ideas  are  repeated.  Consider  a 
particle  with  mass  number  A  and  atomic  number  Z  with  energy  Zecp.  The 
velocity  of  this  particle  is  then  v  =  V 2  Ze<p/Am^  where  m^  is  the 
proton  mass.  The  magnetic  field  B  required  to  turn  the  particle  in  a 

O 

circle  of  radius  r  as  given  by  the  Lorentz  force  is  A  m^  v  /r  =  (Ze/c) 

2 

X  vB.  There  is  then  a  linear  relationship  between  B  and  e<p.  The 

2  2  2 

slope  is  given  by  (2  rn^c  /e  r  )  X  (A/Z).  In  the  mass  spectrometers 
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used  here,  r  was  kept  constant  and  ions,  assumed  monoenergetic,  were 
detected  after  being  bent  through  l80°  as  B  was  varied.  Flux  peaks 
were  observed  at  certain  values  of  B  for  given  cp  and  A.  By  plotting 
these  critical  B  values  against  «p  ,  the  slope  of  the  lines  were  deter¬ 
mined.  Then  the  ratio  of  the  slopes  is  equal  to  the  ratios  of 
(A/Z)  /(A/Z).  where  i  and  j  refer  to  specific  ion  species.  The  slopes 
allow  identification  of  the  ion  species,  but  no  information  about  the 
density.  The  density  comes  from  the  area  under  the  peak  observed 
with  the  spectrograph.  This  area  was  approximated  by  using  the  full 
width  at  half  maximum  (FWHM)  multiplied  by  the  maximum  amplitude  of 
the  peak. 

The  other  method  involved  the  use  of  particle  bursts  called 

21 

pseudo-waves.  Pseudo-waves  were  first  seen  by  Alexeff  et  al.  in 
1968.  These  waves  were  launched  as  a  result  of  negative  potential 
pulses  applied  to  a  grid.  They  explained  the  launching  of  pseudo - 
waves  from  a  grid  by  considering  an  ion  falling  through  a  negative 
potential  well,  produced  by  the  grid  wire.  If  the  pulse  was  removed 
while  the  ion  was  at  the  bottom  of  the  negative  well,  then  the  ion 
will  travel  on  with  a  velocity  v  =  (2  e?  /Am  )£  where  cp  is 
the  maximum  amplitude  of  the  pulse  and  a  singly  charged  ion  is 
assumed.  Likewise,  if  the  particle  is  somewhere  else  in  the  well  the 
velocity  is  v  =  J 2  ecp/Arn^  where  cp  is  the  local  pulse  amplitude.  Also, 
positive  pulses  will  work  as  a  pusher  of  ions  away  from  the  grid. 

The  reason  ion  species  identification  is  possible  is  that 
heavier  ions  will  travel  slower  so  that  the  ion  species  will  become 
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separate  far  from  the  launcher.  This  characteristic  is  inherent  in 

all  waves,  however,  only  pseudo-waves  will  shew  a  change  of  velocity 

with  pulse  amplitude.  If  you  look  at  the  current  collected  by  the 

receiver  the  time  interval  between  the  direct  coupled  signal  and  the 

pseudo-wave  will  decrease  as  <p  is  increased  and  vice  versa.  So 

max 

the  identification  of  pseudo -waves  is  possible  and  the  identification 
of  a  pseudo -wave  for  a  particular  species  is  possible.  The  density 
is  calculated  the  same  as  for  the  mass  spectrometer. 

Experimental  Results 
External  Mass  Spec. 

The  external  spectrometer  was  approximately  nine  inches  from 
the  plasma  anode,  and  the  beam  line  connecting  the  spectrometer  to 
the  chamber  was  sheathed  in  n -metal  to  eliminate  the  fringing  field 
in  the  tube.  The  gyroradius  was  4  cm,  and  the  maximum  magnetic  field 
strength  was  in  excess  of  2  kilogauss.  The  electromagnet  was 
powered  by  a  programmable  power  supply,  which  was  swept  by  a  sawtooth 
from  an  oscilloscope  (Fig.  l8) . 

The  spectrographs  were  recorded  on  an  x-y  plotter  where  the 
x  axis  was  the  output  of  the  programmable  power  supply,  and  the  y 
axis  was  the  output  of  an  electrometer  which  showed  the  ion  current 
collected. 

A  typical  spectrograph  is  shown  in  Fig.  (19)  and  the  magnetic 
field  for  each  ion  species  is  marked  on  the  respective  x  axis.  These 
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values  for  B  are  used  in  the  B  versus  c p  plots  to  identify  the  ion 
species  (Fig.  20). 

Knowing  A/Z,  and  assuming  Z  =  1  for  a  hydrogen  plasma,  allowed 
+  +  + 

identification  of  ,  and  .  The  quantity  of  interest  is 

+  +  ,  V 

the  ratio  of  fluxes  for  and  .  The  ratio  (R)  is  found  by 
measuring  the  areas  under  the  and  peaks  (Fig.  19).  From 
several  spectrographs  it  is  possible  to  see  the  variation  of  R  with 
<PE  (Fig.  20). 

Depending  on  the  choice  of  pQ  and  cp  the  graph  of  R  versus 

<PE  could  be  made  to  do  anything.  However,  taking  cp  -  cp^  =  90  V 

-5 

and  pQ  =  U-  X  10  Torr,  the  curve  of  the  ratio  follows  very  well 
with  the  curve  for  plasma  electron  density  (Fig.  7). 

There  are  several  difficulties  in  using  an  external  spectrom¬ 
eter.  The  beam  line  from  the  chamber  to  the  spectrometer  places  the 
exit  port  in  the  middle  of  a  magnetic  cusp  line.  The  influence  of 
the  surface  magnetic  field  could  greatly  alter  the  ion  extraction 
at  low  ion  velocities  and  would  be  different  for  different  masses. 
Also,  the  long  path  length  allows  the  possibility  of  charge  exchange 
collisions  and  a  change  in  the  ion  spectra. 

Because  of  these  possible  errors  a  method  for  determining  the 
ion  species  composition  needed  to  be  carried  out  inside  the  chamber. 

Pseudo -Waves 

The  pseudo-waves  were  launched  from  a  grid  with  a  diameter  of 
1.8  cm  with  60  lines/inch.  The  collected  current  was  displayed 
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on  an  oscilloscope  using  a  high  gain  plug-in.  The  pseudo-waves  were 
located  hy  varying  the  pulse  amplitude  and  locating  a  wave  that 
moved  properly.  Once  the  waves  were  found,  photographs  were  taken 
of  the  display  screen  for  measurements. 

During  the  analysis  only  one  pseudo-wave  was  observed.  This 
pseudo-wave  had  a  velocity  that  an  H^+  ion  would  have  if  launched 
from  the  tail  of  the  pulse  sheath.  A  pseudo-wave  for  a  H^+  ion  was 
not  observed  because  the  velocity  of  such  an  ion  would  be  buried  in 
the  directly  coupled  signal  of  the  pulse,  or  if  the  distance  traveled 
by  the  wave  was  such  that  the  wave  was  not  buried  then  damping  would 
kill  the  wave  before  it  could  be  detected.  Due  to  the  pseudo-wave 
not  being  observed  an  internal  mass  spectrometer  was  built. 

Internal  Spectrometer 

The  internal  spectrometer  constructed  was  based  on  the  design 
22 

of  Ehlers  et  al.  This  spectrometer  could  produce  a  magnetic  field 
in  excess  of  1  kilo-gauss  with  a  fringing  field  below  100  Gauss.  The 
field  was  produced  by  two  coils  with  each  coil  made  of  550  turns  of 
28  gauge  enamel  insulated  wire.  Between  the  pole  faces  a  copper  box 
was  used  to  short  out  electric  fields  where  the  ion  current  was 
measured.  The  internal  spectrometer  had  a  1.5  cm  gyroradius  and  a 
two  grid  entrance  aperture.  The  outer  grid  was  used  to  keep  out 
primary  and  plasma  electrons,  the  inner  grid  was  used  to  retard  the 
ion  energy.  The  ion  energy  had  to  be  reduced  so  the  magnetic  field, 
lowered  by  heating,  could  bend  the  ion  orbit  in  the  required  semi¬ 


circle 


This  measurement  was  also  inconclusive  due  to  the  poor  resolu 
tion  of  the  spectrometer.  The  poor  resolution  was  a  result  of  the 
fringing  field,  and  that  the  particle  orbit  was  poorly  defined 
because  the  entrance  aperture  was  poorly  collimated. 

The  measurement  of  the  fractional  increase  in  H^+  could  not 
be  completed.  The  methods  are  theoretically  valid  but  technically 
hard  to  use.  The  measurement  is  important  to  show  that  the  observed 
increases  in  the  ion  density  are  due  to  increased  confinement. 


VII .  SUMMARY 


This  thesis  has  presented  experimental  results  of  electro¬ 
static  plugging  of  multidipole  cusps.  The  results  show  increased 
ion  and  electron  density  and  increasing  plasma  potential  with 
increasing  electrode  bias .  A  model  is  presented  and  shown  to  give 
good  qualitative  agreement  with  the  experiment.  This  work  has  also 
found  and  measured  the  contribution  of  secondary  electrons  to  the 
plasma. 

This  work  allows  other  experimentalists  the  ability  to  have 
"high"  density  plasmas  at  elevated  plasma  potentials.  This  freedom 
allows  further  investigation  of  secondary  electrons,  increasing  the 
density  of  a  particular  ion  species  over  another,  or  any  other 
project  requiring  this  combination  of  parameters. 
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Figure  1.  Electrode  and  magnetic  field  conf iguration  in  the  e-pot  devide.  Filament  and 

electrode  biases  voltages  are  shewn  schematically.  Magnetic  poles  face  the  plasma 
and  alternate  in  polarity  as  indicated.  The  magnet  rows  and  electrodes  are  con¬ 
tinued  across  the  top  and  bottom. 
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Figure  3.  Typical  Langmuir  probe  characteristic  (without  secondary  electrons).  The  lines 
A  and  8  subtract  different  amounts  of  primary  electron  current. 
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Figure  4 
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.  Semi -log  plot  of  probe  current  vs.  probe  voltage  (from 
Fig.  3).  The  effect  of  incorrectly  removing  the  primary 
electrons  is  seen  to  raise  the  plasma  electron  tempera¬ 
ture. 


<pE  (VOLTS) 


O 


Figure  5 
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Plasma  potential  is  the  interior  of  the  device  versus 
electrode  bias  at  L  x  10  ^  Torr  with  12  mm  electrode 
separation. 


Figure  7.  Plasma  density  versus  electrode  bias  at  several  values 

of  neutral  pressure  for  6  mm  and  12  mm  electrode  separa¬ 
tion  (same  data  as  Fig.  6). 
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Figure  10.  Plasma  potential  and.  density  versus  electrode  bias 
<%  at  1  X  10  ^  Torr  in  the  presence  of  a  large 
extraction  grid.  Data  is  given  for  floating  and 
grounded  grids. 
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Figure  15.  Graph  of  Eq.  (26). 
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of  Equation  (46). 
of  Equation  ( 47 ) . 
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Figure  15.  Grid  floating  potential  versus  electrode  bias  at  several 
neutral  pressures. 
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Figure  l6. 


A  representative  Langmuir  probe  trace  (labelled  DATA). 

The  plasma  potential  was  81.5  V.  Data  has  been  displaced 
slightly  to  allow  comparison  with  the  fit  to  the  data. 

The  contributions  from  the  plasma  electrons  Ig,  primary 


electrons  I 

P 


indicated. 


and  secondary  electrons  I  are  separately 

s  -6 

The  neutral  pressure  was  5  X  10  Torr. 
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MASS  SPECTROMETER  CIRCUIT 


Figure  19-  Schematic  of  externally  mounted  mass  spectrometer 
experiment. 
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Figure  20.  Representative  spectrograph  of  a  hydrogen  plasma 
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Figure  21.  Top  three  curves  show  the  square  of  magnetic  field  for 
maximum  flux  vs.  plasma  potential.  The  bottcm  graph 
shows  variation  of  density  of  atomic  hydrogen  to 
molecular  hydrogen  with  plasma  potential  when  E  =90' 
and  p  ~  k  X  10  '  Torr.  P 
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